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J> ' Abstract 

It has long been established that the isospin symmetry is slightly broken in the 
nuclear interaction. On the quark level, most studies on the violation of isospin 
f^*) ' symmetry focus on the mass difference between u and d quarks, in addition to 

the (direct) Coulomb interactions between quarks. However, it has been demon- 
strated that there is other source of isospin violation on the quark level. Maltman, 
Stephenson, and Goldman(MSG) showed that the interference effects of QED and 
QCD gave significant contributions to the baryon isomultiplet mass splittings and 
binding energy difference between 3 H and 3 He. In this work, we present the results 
of the effects of these new mechanisms on isospin violation in NN 1 So scatterings 
with the new sets of strength parameters obtained by fitting to the mass splittings 
in the baryon isomultiplets listed in the most recent PDG compilation. We also give 
the matrix elements of the various potential operators which would be useful in the 



study of these charge dependent effects in finite nuclei. 



1. Introduction 

It is well established that isospin symmetry is slightly broken in the nuclear interac- 
tion. However, isospin violation still remains as one of the lesser understood aspects of the 
nuclear force. Both the charge independence breaking (CIB) and charge symmetry break- 
ing (CSB) effects have been extensively studied within the hadronic picture (for recent 
reviews, see Ref. [[[]). For example, n and tt° mass difference in the one-pion-exchange 
potential (OPE) and the exchange of a 771 pair JJ, |3], |J break the charge independence, 
while meson mixings like n — rj and p — u induce charge symmetry breaking interactions 
i 1,0, &!• 

In addition to the interactions induced by the exchange of heavier mesons, the short- 
range isospin violating force could also arise from the quark effects. Well-known sources 
of isospin violation on quark level include the u — d quarks mass difference and the electro- 
magnetic interaction (one-photon-exchange) between quarks. Their effects on the charge 
dependence of the NN scattering length were investigated within the constituent quark 



model in [10, 11, 12] and found to be important but model dependent. The u — d mass dif- 



ference effects on Okamoto-Nolen-Schiffer anomaly |T3[ and isovector mass shifts in nuclei 



14| have also been studied and found to be large. Several years ago, Maltman, Stephen- 
son, and Goldman (MSG) Jl5| proposed a new source of isospin violation on the quark 
level, namely, the interference effect of QED and QCD, i.e., the penguin-like and box-like 
diagrams as depicted in Fig. 1 and 2. They demonstrated that this mixed QED-QCD 
effect gave significant contributions to the baryon isomultiplet mass splittings |TJ| and 



binding energy difference between 3 He and 3 H fL6| |. Another isospin violation interaction 



on the quark level is associated with the quark-pion coupling. Besides the conventional 
long range one-pion-exchange potential (OPEP), the quark-pion coupling also generates 
short-range repulsion and medium-range attraction from the exchange of quarks from 



two different nucleons ||17|| . When the mass difference between charged and neutral pions 
Am n = m n ± — m n o is taken into account, it induces a charge dependent quark-exchange 
OPE (CDEOPE) interaction. 

The difference between proton-proton and proton-neutron singlet S-wave ( 1 5'o) scatter- 
ing lengths is a good testing ground to study isospin violation because the 1 S'o two-nucleon 
system is nearly bound so that the scattering length is large and very sensitive to the de- 
tails of the NN potential. We have studied contributions of the QED-QCD mixing effect 
and CDEOPE interaction to the charge dependence of the scattering length difference in 
1 5'o NN scattering within a nonrelativistic constituent quark model and the preliminary 



results were reported in Refs. [|1| |19 |. In this paper, we present the details of our cal- 
culation and the results of a more extensive study on these quark effects on the charge 
dependence of the scattering length difference in 1 Sq NN scattering. Furthermore, the 
new sets of strength parameters of the QED-QCD mixing potential, obtained by fitting to 
the mass splittings in the baryon isomultiplets listed in the most recent PDG compilation 



2D| , are used in the current calculation. We also give the matrix elements of the various 



potential operators which would be useful in the study of these CD effects in finite nuclei. 

In Sec. 2, we discuss the quark mediated charge dependent interaction. We focus on 
QED-QCD mixing effects and the charge dependent potential associated with the pion 
exchange between quarks when the ^ — n° mass difference is taken into account. In 
Sec. 3 the determination of parameters in the QED-QCD mixing effects is discussed in 
details and the new sets of strength parameters are obtained with the use of a new mass 
splitting value of £+ + — S°. The resonanting group scattering equation used to estimate 
their effect on the AW S-wave( 1 S'o) scattering length is derived in Sec. 4. The results are 
presented and compared with the experimental values in Sec. 5 and the summary is also 
given there. 

2. Quark mediated charge dependent AW interactions 

The aim of this paper is to investigate the various quark effects on the isospin viola- 
tion in 1 Sq AW scattering. Quark-mediated isospin violation mechanisms considered here 
include u — d quark mass difference, quark-exchange Coulomb interaction, QED-QCD 
mixing effects and pion mass difference in quark-exchange OPE interaction. Each effect 
will be individually discussed below. 

2.1. Up-down quark mass difference effect 

The mass difference Sm q = m^ — m u between u and d quark induces isospin violation 
of AW interaction because it explicitly breaks the SU(2) isospin symmetry of QCD La- 
grangian. In the constituent quark model, quark mass difference not only breaks isospin 
symmetry through the kinetic energy of constituent quarks, but also through one-gluon- 
exchange (OGE) interaction between quarks. In 1 5'o, only the Fermi contact hyperfine 
interaction contributes to the isospin violation of AW scattering: 

V 9 m = E -(- ■ ^)(^ ■ ^)-^<s (3) (%), (1) 

iijLi 2 2 A 2 2^3771^- v 3h v ; 

where a s is the strong coupling constant. The AA" splitting m-A — rriN = „ 4 ^ 2 — 260 
MeV leads to a s =1.624, where mo = 330MeV is the average mass of u and d quark. 

2.2. Exchanged pion mass difference effect 

Another important interaction between the quarks arises from the quark-pion coupling. 
Quark-pion coupling is known to be very important to ensure chiral symmetry in quark 
models. It has been used to study various aspects of the AW interaction, for example, 
the conventional long range one-pion-exchange potential (OPEP). In addition, it can also 
generate short-range repulsion and medium-range attraction from the exchange of quarks 
from two different nucleons [0. When the mass difference between charged and neutral 



pions Am n = m n ± — m^o is taken into account, it gives rise to the following charge 
dependent one-pion-exchange interaction between quarks i and j, 

If 2 3 3 e -rn c r 

± i J irqq 



VST(Am,) = - l !g )M . m g ( . 9j) + ll + — + — m — 

3 3 ,„ ,e" mcr 



- T iz T jz {m c {(Oi ■ (Tj) + [H 1 — jSy}- 

3 3 g-mor 

- m {(^-^) + [l + + -2-2 ]Sy} }}, (2) 



where m c = m n ±, itiq = m^o and Sy = 3(c?j • %)(<?« • fy) — c?j • Oj and r = |fy|. 
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Identifying the folded OPEP with the familiar OPEP between nucleons gives f% 



m 



_9_ f2 
25- 



ft NN e-*2 0, with fl NN /ln = 0.079. 



2.3. Exchanged electromagnetic quark effect 

The electromagnetic interaction is an important mechanism of isospin violation. The 
one-photon-exchange diagram induces both Coulomb potential and hyperfine interaction: 



V r = Vf + V*» f , (3) 



where 



V c = J2 a ™Q®^ ( 4 ) 

V- = E-ftftlffl^lr,). (5) 

Here Q u = 2/3, Qd = —1/3 and a em = 1/137. Note that direct Coulomb effect should be 
subtracted. However Coulomb effect between two exchanged quarks has to be taken into 
account. 

2.4. QED-QCD mixing effect 

The QED and QCD mixing effect, as suggested by MSG, refers to those of the penguin- 
like and box-like diagrams as depicted in Figs. 1 and 2, respectively. The penguin-like 
diagrams are the electromagnetic vertex correction to the one-gluon exchange graphs. 
The box-like diagrams are attributed to those of the gluonically dressed versions of the 
basic one-photon exchange graphs. In general, both of them include graphs to all orders 
of a s while only the lowest-order graphs in («a s ) are shown in Figs. 1 and 2. MSG 
based their argument on two chiral constraints on the pion electromagnetic (EM) self- 
energy and demonstrated that these additional terms should be present. The first chiral 



constraint is that [[H|] the contribution from up-down quark mass difference, 8m q , to the 
7r =l= — 7T° splitting is of order 0((Sm q ) 2 ) and small (less than 0.17 MeV) for physical light- 
quark masses. The 7r =l= — it splitting is, therefore, due essentially entirely to pion EM 
self-energy The second is, in the chiral limit |22| ((5m^) em = 0, and the correction to 
the chiral expansion of this constraint is (5m 7T o) em < 0.5 MeV ||23| . Hence the entire 
tp± _ Ti-0 splitting should come from 7T ± EM self-energy However, this cannot be the 
case within the constituent quark model if the one-photon-exchange potential is taken 
to be the only mechanism of electromagnetic interaction between quarks. The opera- 
tor which corresponds to one-photon-exchange graph has the form of Q1Q2F, where F 
includes all the other color, spin and spatial dependences. The expectation values of 
the charge dependent operator Q1Q2F in 7r =l= and n°, are respectively, 2x/9 and — 5x/9, 
where x = {F) n . The chiral constraint of (5ml. ) em = would lead to (5m n ±) em = and 
(m„.± — m n o) < 0.17 MeV, in contradiction with the experiment. On the other hand, if 
we insist in (m ff ± - m^) em ~ (m„± - m^o) exp , then (5m„o) em ~ -f (m^± - m^) exp ^ 0, 
and is inconsistent with the chiral constraint. To resolve this problem, MSG asserted that 
it was necessary to include the penguin-like and box-like diagrams. The explicit forms of 
their effective potentials and associated coefficients are discussed in Sec. 3. 

3. The Hamiltonian 

Our calculation is performed in a non-relativistic constituent quark model with a six- 
quark Hamiltonian, constructed as a sum of one-body kinetic energy, mass and two-body 
potential energy [ID|. It is given as: 



H = K + M + V conf + V glu + V™ + V, + V box + V pengum . 

Since the constituent quark is heavy, we use the non-relativistic form for the free energy 
of the quarks and write 

K + M = ■&£.)+«* 

1=1 l 

Here, rrii is the mass of i-th quark and we use m u = 330MeV. The confining potential is 
taken to be 

Vconf = E -(y4 )( I^' + eo) ' 

i<j=l 

For simplicity, the anharmonic component is not included in V con f. The oscillator length 
a = (3m K) 1/ ' i = 320MeV = 0.617/m" 1 , and the constant e = -385MeV are all fixed 
by the hadron spectra. 

The other terms V g i u , V q °^ ep ', V 7 , are the charge dependent two-body potential terms. 
Their forms and parameters have already been given in Sec 2. Following we only discuss 



the terms associatedwith QCD-QED interference effects. Four effective potential opera- 
tors are introduced by MSG to account for the effects of those penguin-like and box-like 
diagrams: 

/ box 
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2. ..8) are the color SU(3) group generators with standard normalization, 
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A, B, C, D are the strength parameters determined [15] from fit 



Tr{XiXj) 

ting to the isomultiplet mass splittings in the baryon octet and decuplet, and £+ + — E° 
with the assumption that the splittings are induced by the u — d quark mass difference, 
one-photon-exchange potential and the QED-QCD mixing effect discussed above. 



The mass splittings in the baryon octet and decuplet depend only on 8m q , a, [3 
and 7, but not on 5, where a 
(e = (^)nucleon = 1.86MeV, \i = ( 
is necessary to go beyond SU(3)/ baryons, i.e., to systems with different symmetry prop- 
erties and £+ + — S° was chosen in Ref. JT5I for this purpose. The QED-QCD interference 



Ae-Ce 

8TTaem 



3m,m, 



5(3) ( r 



P = Bfi - Dpi, 7 = 4Ce + Dfi, 5 = D/i/A 
ij))nucieon = 1.17 MeV). To determine 5, it 



effects on the charge dependence in NN 1 Sq scattering length reported in [T3|, [H| were 
obtained with A, B, C, D determined in [|TJ with the value of £++ -E° = 0.2 ±0.5 MeV. 
The most recent PDG |2(J compilation now gives S+ + — S° = 0.66 ± 0.28 MeV. We have 
carried out a refitting procedure with this latest value of S+ + — S°. The results are re- 
ported in Sec 5. 



4. The resonanting group method 

The variational scattering equation for the AW system is derived by the resonanting 
group method (RGM). We briefly sketch the method here and refer the readers to Ref. 
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I(J for details. It is based on the Hill-Griffin- Wheeler variational principle, namely, the 



minimization of the functional (ty\H — E\ty), 

(£|\I/> = *(0 = J2 I dRg a (R)® a (£; R), (12) 

a J 

where the 6g state function ^ is a functional of the trial weight function g a (R) which 
depends on the collective variable R, and $ Q (£; R) denotes the 6q wave function in config- 
uration of two clusters A and B. a is a label for the quantum numbers of the clusters and 
£ represents the internal coordinate sets. In this study, we will restrict our consideration 
to the case where both clusters A and B are in the nucleon ground state and write 

$(£; R) = A[<p A <p B ^\X - R)], (13) 

where is the ground state wave function of the 3q cluster and A is the antisymmetrization 
operator, i.e., A = \/^j(l — HitAj&B Pij) i with P^ the permutation operator for particle 
labels % and j. We treat all the interaction terms internal to the clusters, except the 
harmonic potential in V con f, perturbatively up to the first order. This allows us to work 
with the cluster wave functions in the form 

A = 0>-,A)0 c <r, 

3/2 3/2 

m*) = ?*^e-iW+**\ (14) 

where p and A are the conventional Jacobi coordinates in the 3g cluster A and a p = 
(?>mK) l / 4: , a\ = [2>m\K) 1 ^, m\ = 3mm' /(2m + m'), with m the mass of the like quark 
and m! the mass of the unlike quark. <p A , the spatial part of 4>Ai is the ground state 
harmonic-oscillator wave function, <j) c the 3q color singlet sate function and <p aT the flavor- 
spin state function of the nucleon. For simplicity, we will use the charge-symmetric limit 
form of Eq. ([14]) by setting a p = a\ = (3m K) 1//4 . 

The variational principle then gives 

f dR'[H(R,R') - EA{R,R!)\g{R!) = 0, (15) 

where 

H(R,R') = (®(R)\H\<S>(R')) 

= (^(R)\H\^(R')}-9(^(R)\HP lm \^(R')}, (16) 



and 



A{R,K) = ($(R)MR>)) 

= ($ (R)\$%R')}-9($ (R)\P lm \$°(R')}, (17) 



with $°(£;.R) = (J)o(pa,^b)4 > o(pb,^b)8 ( - 3 \X — R). The subscripts l,m of Pi m in Eqs. 
(|16|) and (|17]) are understood to belong to different clusters. Both kernels H (R, R') and 
A(R, R') are divided into a direct (without Pi m ) and an exchange (with Pi m ) part. If we 
define 

(&{R)\HP lm \^{R>)) = H E (R,R>), (18) 

(<S> (R)\P lm \$ (R>)) = K(R,R>), (19) 

then we can express Eq. ([H|) in the form of a Schroedinger equation with a non-local 
potential, after the removal of the CM motion of the two clusters, 



(E-M A - M B ) + V° ep (R) + V°(R) + V°{R)]g{R) + J S(R, R>)g(R)dR> = 0, 



where 

S(R, R!) = -9H E (R, R>) + 9EK(R, R'). (20) 

Ma is the nucleon mass predicted in this model, and fi AB is the reduced mass of A^A^ 
system, 

(o) _ M A M B 1 

with M° = 2m u + m d , M® = 2m d + m u . We have included in the direct term of Eq. (|20|) 
a contribution V® which arises from a— quark coupling. This is to take into account the 
observed medium-range attraction in the A^A^ interaction which has been associated with 
the exchange of low-mass weakly correlated 2ir states. The expressions for all the direct 
terms can be found in fllOl . 



The next main task is to evaluate the exchange kernels of Eqs. (|16D and (p~7|) . Again, 
the expression for K(R, R') and the one-body part of He(R, R') have been presented in 
pTOJI . For the two-body part of He(R,R'), we have here the exchange Coulomb (EC), 



QED-QCD mixing effects and CDEOPE interaction which were not considered in Ref. 
[TTi|| . As in jnj , we split the two-body exchange matrix elements into four components, 



/ , TimPij = F A + T B + T c + Tc' 



l<m 



[Ti.Pi,] + [4(IV + rV,-)^] + [^i'j'Pij] + [(Ti>i» + !>,•„)], (22) 



where the four components, enclosed within brackets, correspond to the four distinct cases 
for the relation between indices /, m of the interacting pair and i, j of the transposed pair: 
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like indices (A), partially like (B), unlike intercluster (C) and the unlike intracluster (C). 
Since the 6g wave function features a simple factorization in color, space and isospin-spin, 
it allows us to discuss independently the overlap matrix element associated with each of 
these variables. Then the only new operators which are not contained in V con f and V g i u 
and not considered before in |Kj are the spatial operators l/Vy and e~ )irij jr^. Only the 
exchange matrix elements of e~^ Tii jr^ are presented in Table 1 since 1/r^ is a special 
case of e~ Mry /ry by letting jjl = 0. We can now write 



S(R,R') = 9[(E - M» A - M B )K(R,R) - !C*(R,R) -V* nf (R,R) -V* U (R,R' 
-Vy (-R, R) - V qqOPEP (R, R') - V box (R, R') - V penguin (R, R% 

In the partial wave expansion, Eq. ([15| ) becomes 



[ 2^% K dR? R 2 

-(E - M A - M B )\ gi (R) + / dR'RR'Si(R, R')gi(R') = 

Jo 

The above integro-differential equation is solved iteratively. At large distance, gi(R) 
should approach 

gi(R) R -^ (kR)(cosSi(k)ji(kR) + sinSi(k)m(kR)) 

~ sin(kR--h + Si(k)), (24) 

where k 2 = 2/i AB (E — Ma — Mb) defines the relative cm momentum. The effective-range 
parameters are calculated from: 

k ■ cot S (k) = -- + -r k 2 + C(A; 4 ), 
a 2 



where a is the scattering length and r the effective range. To obtain the values of Sa, we 

lexp 5 ( 1 / a ) With a exp = a(pp) exp 



use the relation 5a = — a 2 '5(1/ a) with a exp = a(pp) exp ~— 17.0/m. 



5. Results and Discussions 

The results of 5a>csB = o, pp — a nn and SaoiB = \( a pp + a nn) — a np are summarized 
in Table 2-6. Here we have carried out a refitting procedure with this latest value of 
E+ + — S° = 0.66 ± 0.28 MeV given by he most recent PDG |20| compilation. Comparing 
to the previous results |L8L [Tj| where E+ + — E° = 0.2 ± 0.5 MeV was used, it is found 



that SacsB decreases about 0.27 fm and SaciB decreases about 0.15/m. It indicates that 
QED-QCD interference effect exhibts only mild dependence on the value of Ej!" + — E°. 



The range of parameters corresponds to —4.40 MeV > AS*(= E* + — £* ) > 
-5.8 MeV and 5.51 MeV > 5m q > -2.97MeV. On the quark level, the values of 
total SacsB, ranged from 1.35/m to 1.56/ra, are relatively insensitive to the choice of 
parameter sets. On the other hand, the values of SaciB are much more sensitive to the 
choice of parameters. To understand this difference one needs to study each isospin vio- 
lation mechanism in more details. 

The dominant CSB effects on the quark level are QED-QCD interference and quark 
mass difference effects. Although individually each effect is very sensitive to the chosen 
parameters, (Note that quark mass difference effect is ranged from — 1.83/m to 3.24/ra 
with respect to different parameter sets; the penguin-like diagrams contribution is ranged 
from 3.91/ra to — 1.32/m and the box-like diagrams contribution is ranged from — 1.14/ra 
to — 0.74/m) their sums still keep approximately constant. 

On the other hand, similar cancellations do not occur in the CIB case. Since V(Sm q ) 
and penguin-like diagrams are charge independent, the sensitivity of 5clcib with repsect to 
the chosen parameters is completely due to the fact that the box-like diagrams contribute 
to SaciB ranging from —2.82fm to — 1.82 fm. Furthermore, the CIB mechanisms on the 
quark level which are independent of chosen parameter sets include the quark-exchange 
Coulomb interaction and quark-exchange OPE interaction. The quark-exchange Coulomb 
(EC) effect is about 0.38/ra, the quark-exchange OPE effect (CDEOPE) is only 0.08/ra. 
Their sum is only about 0.46 /ra, much less than the box-like diagrams contributions. It 
explains why our values of SaciB are more sensitive to the parameters. 

Furthermore, let us compare our results with the experimental data. The 1 Sq proton- 
proton scattering length is very accurately measured, but the subtraction of the direct EM 
interaction is model dependent J24|. A commonly quoted value is a pp = —17.0 ± 0.2 fm. 
The corresponding value for neutron-neutron scattering as determined from n~d — ► 'ynn 



is a nn = —18.5 ± OAfm [25]. Our values of quark contributions to the CSB range from 
1.35/m to 1.56/ra, are very close the experiemtnal value. However, Coon and Niskanen 
|26|j found that from CSB vertex corrections and mass differences of the intermediate 



baryons, the two-pion exchange would contribute to SacsB- This effect would be as large 
as 1.37 fm. The sum of the meson and quark contributions would then be larger than the 
experiemtnal value of 1.5/ra, but is still within the uncertainty incurred in the subtraction 



of direct EM interaction in pp scattering [E| 



The 1 5'o neutron-proton scattering length is — 23.78 ±0.001/ra, therefore a(pp — np) = 
8a>ciB + \8acsB — 6.78 fm |27f . It was claimed in Ref. J28| this result can be explained 



by various meson-mediate mechanisms. The dominant contribution is from the pion mass 
difference effect in OPEP, its size was 3.6 ± 0.15/m. Another effect is due to the jn 
exchange which was 1.1 ± OAfm. Such an optimstic picture can be called into ques- 
tions with a closer look. First, the pion mass difference effect in OPEP only contributes 
to a(pp — np) 2.64 ± 0.16/ra according to Refs. [^|, |30| and confirmed by the calcula- 
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tion of the Ref. |2"9fl . For the 771"- exchange effect mediated by the seagull vertex, Ref. 
@, § obtained only a contribution of 0.67 ± 0.02/m to a(pp — np) through the solution of 
Schrodinger equation with RSC potential. A recent calulation of the same effect even gave 
a result of —0.35 to — 0.53/m with opposite sign [|J. Taking these two together would 
amount to a reductionof about 1.40 to 2.60/m from the total estimate in Ref. [p!S|] . The 
quark contribution to a(pp — np) in our calculation is ranged from —1.68 to — 0.57 fm. It 
clearly makes the agreement between experiments and theories worse since it is opposite 



in sign to the experimental value. If one uses the the estimate of Refs. |3], g9[ |30] for 
the pion mass difference effect in OPEP, Ref Q for 771"- exchange and Ref. [28| for all 
other meson-mediated effects respectively, and adds our values of quark contributions, 
then a(pp — np) ranges from 1.31 to 2.42/m, well below the experimental value. 

Although we have calculated Sclcib (CDEOPE) and found its value only about 0.08/m, 
its contribution may be underesteemed. Stancu et. al. [J3TJ] pointed out that the single 



channel Resonating Group Method maybe not adequate for the operators with this fac- 
tor :(tj • 7jr)((7j ■ Oj). Because the dominant configuration is no longer |s 6 [6] [33]_fs) but 
|s 4 p 2 [42] [51]^,s), the reason is because of its specific flavor-spin symmetry. It requests 
further studies to esteem more accurately on the size of this effect. 

Furthermore, in the previous fit, the 7r-quark coupling inside the hadron has been 
completely neglected since we truncated the 7r-quark coupling in the long and medium 
range. To be consistent, quark-exchange OPE effect should be also included in the fit of 
isomultiples of baryon spectrum if <5ac7g (CDOPE) is taken into account. Such a refit is 
expected to provide completely different values of all parameters here, including a s , 5m q 



and A,B, C, and D. However the recent calculation of Shih |£2] showed the 7r-quark cou- 
plings inside hadron are small effects and refitted 5m q increases only by 0.3MeV, therefore 
we expect the new sets of parameters generated from this refit will not dramatically differs 
with ours. 

In conclusion, the new sets of strength parameters of the QED-QCD mixing potential, 
obtained by fitting to the mass splittings in the baryon isomultiplets listed in the most 
recent PDG compilation |2(J, is used in the current calculation. The effects of interference 



of QED and QCD on the isosp in- violation of 1 So NN scattering are significant, and very 
sensitive to the chosen parameters. The effects of QED-QCD interference should be also 
significant in charge-dependence in nuclear many-body systems such as u and d quark 



mass difference effect does [JL4]], therefore , we give the matrix elements of the various po- 
tential operators which would be useful in the study of these charge dependent effects in 
finite nuclei. The charge dependent NN potential generated from several mechanisms on 
the quark level explains the CSB effect of 1 5'o NN scattering very well, without including 
the questionable p — u mixing effect. But the CIB effect on the quark level carries the 
opposite sign with respect to the experimental value and the agreement between theory 
and experiment is still unavaible, so that the charge independence breaking in the NN 
interaction remains as an open issue in the quark models. 
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Type 


Function 


A 


2 exp(-§/x|R-R'|) 
3 |R-R'| 


B 


||^{exp[-f|3R-R'|]er/( V / |a[||3R-R'|-&]) + exp(f|3R-R'|) 
er/( v /|a[||3R - R| + f^\) - 2sinh(f |3R - R'|)} + (R <-► R!) 


C 


g^{exp(-f|RH-R'|)er/(a[f|RH-R'|-^])H-exp(f|R + R'|) 
er/(a[f |R + R'| + &]) - 2 sinh(f |R + R'|)} + (R <-> R') 


C 


a^-^-erfc(^)exp[^} 



Table 1: The expectation of (^^i*>)/K°(R, R'). K°(R, R')=(§^y/ 2 e xp[-±§a 2 ( J R 2 + 
i?' 2 ) + |a 2 R • R]. To get (j-P% simply take /i=0. 



^fM^V 



vjM^V 



-S^^^X 



>vvw^ 



Aw^ 



>VVW^ 



Figure 1: Penguin-like diagrams for 0(aa s ) isospin violation arising from QED-QCD 
interference. The sharply varying, smoothly varying, and solid lines refer to photons, 
gluons, and quarks, respectively. 



Figure 2: Corresponding box-like graphs. Notations same as Fig. 1. 
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5m q 


box (el) 


box (mag) 


pen (el) 


pen (mag) 


subtotal 


EC 


CDEOPE 


total 


5a C sB 


-1.83 


-0.05 


-1.09 


0.46 


3.48 


0.97 


0.38 





1.35 


SaciB 





-1.63 


-1.19 








-2.82 


0.38 


0.08 


-2.36 



Table 2: The quark contributions (in fm) to SacsB = a PP — a>nn, 8clcib = |(a P p + a nn ) — a np . 
The notation in the entries is, e.g., pen (el) refers to the electric effect induced by the 
penguin-like graphs. EC and CDEOPE refer to the effect of quark-exchange Coulomb and 
OPE interaction. Here AS* = -AAOMeV, 5m q = -2.97MeV, A = -1.96, B = 7.90, 
C = 9.41, D = -13.28. 





5m q 


box (el) 


box (mag) 


pen (el) 


pen (mag) 


subtotal 


EC 


CDEOPE 


total 


5a C sB 


-0.52 


-0.03 


-1.00 


0.32 


2.28 


1.05 


0.38 





1.43 


5a C iB 





-1.48 


-1.09 








-2.57 


0.38 


0.08 


-2.11 



Table 3: The notation is same with Table 2. Here AS* = — A.75MeV , 5m q = —0.85MeV, 
A = -1.78, B = 7.20, C = 6.91, D = -8.71. 
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5m q 


box (el) 


box (mag) 


pen(el) 


pen (mag) 


subtotal 


EC 


CDEOPE 


total 


5a C sB 


0.76 


-0.03 


-0.91 


0.22 


1.10 


1.14 


0.38 





1.52 


5a C iB 





-1.33 


-0.98 








-2.31 


0.38 


0.08 


-1.85 



Table 4: The notation is same with Table 2. Here AS* = -b.lOMeV, 5m q = 1.27 MeV, 
A = -1.61, 5 = 6.51, C = 4.41, D = -4.14. 





5m q 


box (el) 


box (mag) 


pen (el) 


pen (mag) 


subtotal 


EC 


CDEOPE 


total 


Sa C sB 


1.98 


-0.02 


-0.81 


0.10 


-0.10 


1.15 


0.38 





1.53 


5a C iB 





-1.19 


-0.87 








-2.06 


0.38 


0.08 


-1.60 



Table 5: The notation is same with Table 2. Here AS* = —5A5MeV, 5m q = 3.39MeV^, 
A = -1.43, B = 5.81, C = 1.90, D = 0.43. 





5m q 


box (el) 


box (mag) 


pen (el) 


pen (mag) 


subtotal 


EC 


CDEOPE 


total 


Sa C sB 


3.24 


-0.03 


-0.71 


-0.02 


-1.30 


1.18 


0.38 





1.56 


5a C iB 





-1.05 


-0.77 








-1.82 


0.38 


0.08 


-1.36 



Table 6: The notation is same with Table 2. Here AS* = — 5.80Mey, 5m q = 5.51MeV, 
A = -1.26, B = 5.11, C = -0.60, D = 4.99 
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